Galactic turbulence and paleoclimate variability 



A. Bershadskii 
ICAR, P.O.B. 31155, Jerusalem 91000, Israel 

The wavelet regression detrended fluctuations of the reconstructed temperature for the past three 
ice ages: approximately 340000 years (Antarctic ice cores isotopic data), exhibit clear evidences of 
the galactic turbulence modulation up to 2500 years time-scales. The observed strictly Kolmogorov 
turbulence features indicates the Kolmogorov nature of galactic turbulence, and provide explanation 
to random-like fluctuations of the global temperature on the millennial time scales. 

PACS numbers; 98.38.Am, 92.70.Gt, 96.50.Xy, 92.70.Qr 
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The Earth chmate is an open system. Not only the Sun 
is the main source of energy for the climate dynamics but 
also galactic cosmic rays (GCR) can effect the climate 
on global scales (see, for instance, and references 

therein). On the other hand, recent paleoclimate recon- 
structions provide indications of nonlinear properties of 

Earth climate at the late Pleistocene [i3,lil,[il (the 

period from 0.8 Myr to present). Reconstructed air tem- 
perature are known to be strongly fluctuating on millen- 
nial time scales (see, for instance, figures 1 and 2). While 
the nature of the trend is widely discussed (in relation to 
the glaciation cycles) the nature of these strong fluctua- 
tions is still quite obscure. The problem has also a tech- 
nical aspect: detrending is a difficult task for such strong 
fluctuations. In order to solve this problem a wavelet 
regression detrending method was used in present inves- 
tigation. 
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FIG. 2: The wavelet regression detrended fluctuations from 
the data shown in Fig. 1. 
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FIG. 1: The reconstructed air temperature data (dashed line) 
for the period 0-340 kyr. The data were taken from Ref. [ly] 
(see also Ref. [U)- The solid curve (trend) corresponds to a 
wavelet (symmlet) regression of the data. 



Figure 1 shows reconstructed air temperature data 
(dashed line) for the period 0-340 kyr as presented at 
the site 16| (Antarctic ice cores isotopic data, see also 
Ref. [17|). The solid curve (trend) in the figure cor- 
responds to a wavelet (symmlet) regression of the data 
(cf Ref. [i3|)- Figure 2 shows corresponding detrended 
fluctuations, which produce a statistically stationary set 
of data. Most of the regression methods are linear in 
responses. At the nonlinear nonparametric wavelet re- 
gression one chooses a relatively small number of wavelet 
coefficients to represent the underlying regression func- 
tion. A threshold method is used to keep or kill the 
wavelet coefficients. In this case, in particular, the Uni- 
versal (VisuShrink) thresholding rule with a soft thresh- 
olding function was used. At the wavelet regression the 
demands to smoothness of the function being estimated 
are relaxed considerably in comparison to the traditional 
methods. 
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FIG. 3: A small-time-scales part of the autocorrelation func- 
tion defect of the wavelet regression detrended fluctuations 
from the data shown in Fig. 2. The straight line indicates 
the Kolmogorov's '2/3' power law for the structure function 
in In-ln scales. 
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FIG. 4: A small-time-scales part of the structure functions 
Sp (p=2,3) for the wavelet regression detrended fluctuations 
from the data shown in Fig. 2. The insert shows the structure 
functions for p = 0.2, 0.5, 0.7, 1. The straight line indicates 
scaling in the In-ln scales. 



STATISTICAL ANALYSIS 

The wavelet regression detrended fluctuations shown 
in Fig. 2 represent a statistically stationary set of data 
and we can use such standard statistical tools as correla- 
tion function C(t) and structure functions S'p(t) (of dif- 
ferent orders p) for it analysis. The correlation function 
defect 1 — C(t) is proportional to the second order struc- 
ture function 5*2 (t). Therefore, we can compare results 
obtained by these different tools. First let us look at au- 
tocorrelation function C (r) in order to understand what 
happens on the millennial time scales. Figure 3 shows 
a relatively small-times part of the correlation function 
defect. The In-ln scales have been used in this figure in 
order to show a power law (the straight line) for the sec- 
ond order structure function: 52(t) = {\x{t + T) — x{t)\'^) 



1~C{t) cx52(t) oct^/^ 



(1) 



This power law: '2/3', for structure function (by virtue 
of the Taylor hypothesis transforming the time scaling 
into the space one [3],[2^) is known for fully developed 
turbulence as Kolmogorov's power law. 

Although, the scaling interval is short, the value of 
the exponent is rather intriguing. This exponent is 
well known in the theory of fluid (plasma) turbulence 
and corresponds to so-called Kolmogorov's cascade pro- 
cess. This process is very universal for turbulent flu- 
ids and plasmas 2l|. In spite of the fact that magnetic 
field is presumably important for interstellar turbulence 
the Kolmogorov description can be still theoretically ac- 
ceptable even in this area [l^l , [11] , 24 1 . Moreover, the 



to kpc. In order to support the Kolmogorov turbulence as 
a background of the wavelet regression detrended temper- 
ature modulation we calculated also structure functions 
S'p(r) — {\x{t + t) — x{t)\P) with different orders p. In 
the classic Kolmogorov turbulence (at very large values 
of the Reynolds number [19i].[25i]) 
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Kolmogorov-type spectra were observed on the scales up 



(at least for p < 3). Figure 4 shows a small-time- 
scales part of the structure functions Sp with p = 
0.2, 0.5, 0.7, 1, 2, 3 for the wavelet regression detrended 
fluctuations from the data shown in Fig. 2. The straight 
lines are drawn in order to indicate scaling in the In-ln 
scales. Figure 5 shows as circles the scaling exponent 
against p for the scaling shown in Fig. 4 (the expo- 
nents were calculated using slopes of the straight lines - 
best fit, in Fig. 4). The bars show the statistical errors. 
The straight line in Fig. 5 corresponds to the strictly 
Kolmogorov turbulence with — p/3 Eq. (2). One 
can see good agreement with the Kolmogorov turbulence 
modulation. 



DISCUSSION 

For turbulent processes on Earth and in Heliosphere 
the Kolmogorov's scaling with such large time-scales cer- 
tainly cannot exist. Therefore, one should think about a 
Galactic origin of Kolmogorov turbulence with such large 
time-scales (let us recall that diameter of the Galaxy is 
approximately 100,000 light years). This is not surpris- 
ing if we recall possible role of the galactic cosmic rays 
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FIG. 5: The scaling exponent ("p against p for the scahng 
shown in Fig. 4 (circles). The triangles correspond to the 
cosmic microwave radiation (QMASK map) [sil. The straight 
line is drawn in order to indicate the Koimogorov's turbulence 
with Cp =p/3 Eq. (2). 



for Earth climate. Galactic cosmic ray intensity at the 
Earth's orbit is modulated by galactic turbulence pot - 
On the other hand, the galactic cosmic rays can deter- 
mine the amount of cloud cover (a very significant climate 
factor) on global scales through the massive aerosols for- 
mation (see, for instance, Thus, the galactic 
turbulence can modulate the global temperature fluctu- 
ations by the Kolmogorov scaling properties (Figs. 3-5) 
on the millennial time scales. 

If one knows the characteristic velocity scale v for the 
Taylor hypothesis one can estimate outer space-scale 
of corresponding galactic turbulence as L > 2500y x v. 
However, it is not clear what estimate we should take for 
the V. For instance, one could try velocity of the solar 
system relative to the cosmic microwave background 
(CMB) rest frame: v ^ 370 km/sec. In this case one 
obtains L ^ Ipc. It should be noted that in recent 
paper [2^ it is suggested that the typical outer scale 
for spiral arms can be as small as Ipc and in interarm 
regions the outer scale can be larger than lOOpc. Since 
the solar system and Earth are at present time within 
the Orion Arm this suggestion is in agrement with the 
above estimate. Although the suggestion of the Ref. 
[29I is still under active discussion the paleoclimate 
consequences of this suggestion can be very interesting 
and we will discuss one of them here. Namely, when 
orbiting the Galactic center the solar system and 
Earth are in the interarm regions the reverse of the 
Taylor hypothesis provides us with the outer time-scale 
^ 2500 X 100 years. This time scale is larger then 
any known glaciation period (which are determined by 
the periods related to orbiting Earth around Sun, see 
for instance 



to the glaciation cycles to occur when the solar system is 
in the interarm regions. The Earth deglacitaion related 
to the interarm regions was suggested in Refs. 
and explained by difference in intensity of the cosmic 
ray flux in the spiral arms and in the interarm regions. 
It is difficult to estimate at present time which of the 
two mechanisms is more efficient. In any way they are 
working to the same end and both are open to discussion. 

The same turbulence that modulates the gala ctic cos- 
mic rays in the galactic disk and in its halo |24 | . |2 6l can 
also modulate the cosmic microwave radiation [27| . Nat- 
urally, the effect of the galactic modulation of the cos- 
mic microwave radiation is different in different regions of 
the sky and in different frequency bands. An additional 
component of galactic microwave emission from spinning 
dust, for instance, has a peak in intensity between 20 
and 40 GHz The so-called QMASK map of cos- 

mic microwave radiation combines observations made in 
a vicinity of the North Celestial Pole in Ka and Q bands 
(26-36 GHz and 36-46 GHz respectively) The map 
was generated by subdividing the sky into square pixels 
of side 8 ~ 0.3° and consists of 6495 pixels. The data are 
represented using three coordinates x, j/, of a unit vector 
R in the direction of the pixel in the map (in equatorial 
coordinates). Since some pixels are much noisier than 
others and there are noise correlations between pixels, 
Wiener filtered maps are more useful than the raw ones. 
Wiener filtering suppresses the noisiest modes in a map 
and shows the signal that is statistically significant. The 
space increments of the cosmic microwave temperature 
T can be defined as 



AT, = (T(R + r) - r(R)) 



(3) 



where r is dimensionless vector connecting two pixels of 
the map separated by a distance r, and the structure 
functions of order p as (| AT,!^) where (.) means a statis- 
tical average over the map. The scaling 
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(4) 



11|). Strong turbulent fluctuations of the 



cosmic rays flux on such large time-scales should prevent 



was observed for these increments in Ref. [3l| for the 
QMASK map. Corresponding values of are shown in 
Fig. 5 as triangles. 
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